Mach cones in an evolving medium 
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The energy and momentum lost by a hard parton propagating through hot and dense matter has to 
be redistributed in the nuclear medium. Apart from heating the medium, there is the possibility that 
collective modes are excited. We outline a formalism that can be used to track the propagation of 
such a mode through the evolving medium if its dispersion relation is known. Under the assumption 
that a sound wave is created, we track the jet energy loss as a function of spacetime and follow the 
resulting mach cone throughout the fireball evolution. We compare with the angular correlation 
pattern of hard hadrons as obtained by the PHENIX collaboration and find good agreement with 
the data provided that a substantial fraction of jet energy (~ 90%) is deposited into a propagating 
mode and that the hot matter can be characterized by an EOS with a soft point (not necessarily a 
mixed phase). 

PACS numbers: 25.75.-q 



I. INTRODUCTION 



Announcements have recently been made by the four de- 
tector collaborations at RHIC pj that a new state of 
matter, distinct from ordinary hadronic matter has been 
created in ultrarelativistic heavy- ion collisions (URHIC). 
A new and exciting challenge for both experiment and 
theory is now to study its properties. The energy loss of 
hard partons created in the first moments of the collision 
has long been regarded a promising tool for this purpose 
13 E10 00] However, before this tool can be fully 
exploited, the influence of the medium evolution on the 
energy loss has to be understood in a more quantitative 
way (cf. 0). 

Recently, measurements of two-particle correlations in- 
volving one hard trigger particle have shown a surprising 
splitting of the away side peak for all centralities but 
peripheral collisions, qualitatively very different from a 
broadened away side peak observed in p-p or d-Au colli- 
sions [3 . Interpretations in terms of energy lost to prop- 
agating colourless 0,0] and coloured ^2| sound modes 
have been suggested for this phenomenon. In a 
Mach cone structure has been demonstrated to be vis- 
ible in the angular pattern of photon emission if a source 
term for the medium energy is inserted in a h ydr o code, 
although this is at odds with the results of where 
this form of the source term does not yield a propagating 
Shockwave. In addition, we note that energy densities 
more than 20 times above the value of the surrounding 
medium in the vicinity of 1 fm around the jet as shown 
m Fig. 1 of [13 do not seem to reflect the situation re- 
alized at RHIC, a possible reason for the discrepancy in 
the finding between [13 and [13. 

In the following, we investigate if such a mechanism can 
in agreement with the measured hadronic data provided 
that a realistic evolution model of the hot and dense mat- 
ter is used. 



II. SHOCKWAVE EXCITATION, 
PROPAGATION AND FREEZE-OUT 

Throughout this manuscript we focus on central collisions 
only. Since we are interested in the deposition of lost jet 
energy into the medium, our first task is to determine 
the spacetime pattern of energy loss. 
For the description of the evolution, we base our inves- 
tigation on the fireball model outlined in ^3] and [Tfij - 
The model is tuned to describe transverse mass spectra 
of hadrons, HBT correlation radii and hadronic dN/dy 
distributions. The model has successfully predicted the 
photon yield 0, ITU fl8| at RHIC and gives a fair de- 
scription of Raa m the high region without requiring 
unrealistically high gluon multiplicities or transport coef- 
ficients far from the perturbative estimate [3 ■ Therefore 
the model can be expected to give a fair representation 
of the relevant physics in both early and late stages of 
the evolution. 

For the energy loss calculation we use results of ^3 to 
obtain the probability P(AE) to lose the amount of en- 
ergy AE from the two key quantities, plasma frequency 
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and averaged momentum transfer 
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in a static equivalent scenario which are calculated along 
the path of the hard parton through the medium. Here 
we take = ce(£) 3 / 4 with e the local energy density 
of the medium as determined by the fireball evolution 
[l4| and c = 3 as required for the description of Raa [3 
and close to the perturbative c ~ 2 j^j ■ Since we are not 
interested in folding the result with a steeply falling spec- 
trum but rather into the energy deposited on average in 
a given volume element we focus on the average energy 
loss (AE) = J °° P(AE)AEdAE in the following. For 
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illustration of the resulting loss pattern, we show the en- 
ergy deposition for a hard parton at midrapidity (rj = 0), 
in the transverse (x, y)-plane at y = and propagating 
into positive x direction as a function of initial position 
x in Fig. ^ 



and propagate the part of the shockfront remaining in the 
midrapidity slice (i.e. in the detector acceptance). Each 
piece of the front is propagated with the local speed of 
sound and the angle of propagation is continously ad- 
justed as 
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FIG. 1: Energy deposition from hard quark at r\ = 0, y = 
propagating into positive x as function of proper time in the 
CM frame r for different initial positions x. The shape of the 
energy loss pattern for gluons is qualitatively similar but has 
a larger absolute normalization. 

The result shows the initial rise of the differential energy 
loss as an approximately linear increase reminiscent of 
the L 2 dependence of the energy loss in a static medium 
and a rather sudden drop in dE/dx once expansion signif- 
icantly dilutes the medium. Once transverse expansion 
starts becoming important, the high power of the volume 
growth with time wins out over the advantage of having 
a longer decoherence time. It is the finite formation time 
of radiated quanta which implies that lost energy does 
not appear in the medium immediately which causes the 
curves to rise from zero. 

Our further investigation relies now on the assumption 
that a fraction / of the energy lost to the medium excites 
a collective mode of the medium. Since our framework 
makes use of energy-momentum conservation laws, we 
do not have to address the important (though unsolved) 
question by what mechanism such a propagating mode 
is excited. We considered the case where colorless sound 
characterized by a dispersion relation E = c 2 p is excited 
by that energy fraction /. The remaining energy fraction 
(1 — /) in essence heats the medium and leads to some 
amount of collective drift along the jet axis to conserve 
longitudinal momentum. 

We calculate the speed of sound c s locally from a quasi- 
particle description of the equation of state as fitted to 
lattice results as c s = dp{T)/de(T). This EOS shows 
a significant reduction of c s as one approaches the phase 
transition but does not lead to a mixed phase. The dis- 
persion relation along with the energy and momentum 
deposition determines the initial angle of propagation of 
the shock front with the jet axis (the 'Mach angle') as 
<t> — arccos c s . We discretize the time into intervals At, 
calculate the energy deposited in that time interval E(t) 



f TE c.(r)dT 
(r - t e ) 
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where c s (t) is determined by the propagation history. 
Since sound propagates in the local restframe, the shock 
front element is also carried away by the local flow veloc- 
ity. We illustrate this by showing the Mach cone in the 
transverse r\ — plane for a jet travelling into positive x 
direction, originating at x = —6 fm and either y = (i.e. 
going through the fireball center) or y = 3 fm in Fig. [21 
Both the effect of the soft point in the EOS narrowing the 
cone at late times and the distortion of the cone in po- 
sition space are clearly visible. Note that the distortion 
by flow is sizeable in position space since flow velocities 
vt < 0.7 are of the same order of magintude as the speed 
of sound c s ~ 0.3..0.5. However, since the effect of flow is 
already included in the standard calculation of mj- spec- 
tra, the direction of flow of excess momentum contained 
in the Mach cone in momentum space is hardly changed. 
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FIG. 2: Mach cones excited by a jet travelling from x = —5 
fm into positive x direction (with (0,0) the fireball center) for 
y — (left panel) and y = 3 fm (right panel). 

Once an element of the wavefront reaches the freeze-out 
condition T = Tp , a hydrodynamical mode cannot prop- 
agate further. We assume that the energy contained in 
the Shockwave is not used to produce hadrons but rather 
is converted into kinetic energy. In the local restframe, 
we then have a matching condition for the dispersion re- 
lations 



E = c s p and E = yj M 2 + p 2 - M 



(4) 



where M = V (p(Tp) + e(Tp)) is the 'mass' of a volume 
element at freeze-out temperature. 

Once we have calculated the additional boost a volume 
element receives from the Shockwave using the matching 
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conditions, we use the Cooper- Frye formula 
d 3 N 



E- 



= 9 

<Pp (27r) s 



flow shock 

into a hadronic distribution 



danp^exp — =dxS(x,p) 

1 f 

(5) 

to convert the fluid element 
The resulting momentum 
spectrum is thus a thermal two component spectrum re- 
sulting from an integrations involving volume not part of 
the shockwave and volume receiving an additional boost 
from the shockwave. 

From those considerations one can also calculate the typ- 
cial number particles getting an additional boost from the 
shock front. An upper bound for c s at Tp would be the 
speed of sound in a hadronic resonance gas c s f=a \/0. 2 
p2| with a freeze-out temperature of « 110 MeV fl4j . 
From those considerations one finds for an energy in the 
collective mode of about 5 GeV a number of 20 pions 
which are part of the shockwave. Assuming e.g. c s = 0.3 
one finds « 50 pions. The mach cone would also con- 
tain a smaller number of additional kaons and nucleons. 
A 10 GeV jet can therefore potentially involve 0(100) 
particles. 



reproduces the width of the far side peak in the case of 
d-Au and 60-90% peripheral Au-Au collisions. 
With vertex, energy and direction of the away side jet 
fixed, we calculate dE/dr of the outgoing parton. We 
stop the calculation when a significant fraction of the en- 
ergy is lost to the medium. We generate and propagate 
a shockwave by the formalism outlined in the previous 
section and fold the resulting hadron spectra with the 
experimental acceptance cuts. For the remnants of the 
away side jet, thermal physics and local boosts determine 
the spectral excess, therefore we have no ambiguity as- 
soicated with the choice of recombination/fragmentation. 
In each event we assume that a fraction (1 — /) of the 
energy lost from the hard parton heats the medium. We 
assume that due to momentum conservation this will lead 
to additional flow into the direction of the original away 
side parton. Likewise, we account for the possibility of 
a punchthrough if the initial vertex is very peripheral 
and both near and away side parton propagate near- 
tangential to the surface (however we note that these 
events are extremely rare and play almost no role in the 
calculation). 



III. SAMPLING THE MEDIUM EVOLUTION 

In the following we compare to the hard two particle cor- 
relation data. For this purpose we simulate the PHENIX 
trigger conditions as closely as possible using a Monte- 
Carlo approach. We start by generating vertices with a 
distribution weighted by the nuclear overlap 



T AA {h) = J dzp 2 (b,z). 



(6) 



We then determine the jet momentum and parton type by 
randomly sampling partonic transverse momentum spec- 
tra generated by the VNI/BMS parton cascade as de- 
scribed in Calculating the energy loss of the near 
side parton, we decide if the experimental trigger condi- 
tion is fulfilled. Since the experiment triggers on a hard 
hadron in the transition between the recombination and 
fragmentation regime, the model at this point cannot im- 
plement the trigger condition exactly. Instead, we require 
the trigger condition to be fulfilled by the parton and 
have checked that the model results do not change signif- 
icantly when the near side trigger threshold is increased 
by 2 GeV. We note that this procedure places the ver- 
tices fulfilling the trigger condition close to the surface 
of the produced matter, i.e. in our model the medium is 
rather opaque, in agreement with the conclusions of [2j| . 
Once the vertex and momentum of a near side jet has 
passed the trigger condition, we determine the direction 
of the far side parton in the transverse plane. In order 
to take into account intrinsic hr, we do not propagate 
the away side directly opposite to the near side jet but 
allow for a random angle. We have verified that this 
distribution, folded with the width of the near side peak 



IV. RESULTS 

We present the resulting 2-particle correlation structure 
on the away side in Fig. El Since the near side associated 
particle distribution is not calculable in the framework 
outlined above we will only focus on the away side region 
and make no attempt to describe near side correlations, 
hence our results are normalized to the integral of the 
away side correlation structure. We have chosen de- 
grees as the direction opposite to the near side jet. 
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FIG. 3: Away side correlation structure for different scenarios. 

We present three different scenarios: Our best fit of the 
only model parameter / to the data (leading to / w 0.9) 
a scenario where / is set to 0.5 and one in which we do 
not determine c s by the EOS but set it to the ideal QGP 
value c s ~ 1/V3- 

One can read off two conclusions immediately from the 
comparison with the data: First, a large fraction (about 
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90%) of the lost energy has to excite a Shockwave in order 
to reproduce the data. Since most of the 'Mach ring' in 
momentum space is outside the acceptance but any flow 
or punchthrough is predominantly inside, even a moder- 
ate fraction of jet energy heating the plasma would lead 
to a strong signal at zero degrees. In order to find a 
dip, the mechanism transporting energy away from this 
region has to be very efficient. 

Second, the angular position of the peak is capable of 
measuring the speed of sound to some degree. We observe 
that on average it has to be significantly smaller than the 
ideal QGP value and that the data are in agreement with 
the situation as seen on the lattice, i.e. a soft point but no 
mixed phase. However, it would be premature to claim 
that this is the only possibility to describe the data. 

V. SUMMARY AND OUTLOOK 

We have demonstrated that the idea of Shockwaves being 
excited by the energy lost from a hard parton into the 
medium is able to account for the observed splitting of 
the away side peak if a realistic model for the medium 
evolution is used. If this is the underlying mechanism, 
it has to dominate over simple heating of the plasma, 
as this would produce a significant peak in forward di- 
rection. In this scenario, the away side peaks essentially 
represent thermal physics, and we expect the momen- 
tum distribution at the peak between 1 and 2 GeV to 
be dominated by boosted thermal spectra and not by jet 
fragmentation or recombination physics. 
We emphasize that the framework we have used is quite 
general and relies mainly on energy momentum conser- 
vation and the assumption of a deposition of a consider- 
able fraction of the energy in a collective mode. We used 
a specific dispersion relation (colorless sound) to demon- 
strate that the two particle correlations with hard trigger 
have some sensitivity to c s in hot and dense matter. We 
see an average c s far from the ideal QGP value and well 
compatible with a crossover transition with a soft point. 
Recently, measurements of 3-particle correlations have 
been presented by both the PHENIX 24] and the STAR 
PH collaboration. While we plan to address a calculation 



of 3-particle correlations in the framework outlined here 
in a future publication, we'd like to remark at this point 
that since the Mach cone is a collective phenomenon that 
can - as discussed above - involve typically of the order of 
20-50 particles for the given trigger conditions. While all 
of them are strongly correlated with the original away- 
side parton (and hence with the observable near side jet), 
once this correlation is subtracted only the momentum 
balance (with respect to the away-side parton momen- 
tum) Pl = Paway and Pt = remains to account 
for residual correlation inside the cone. With momen- 
tum distributed across a significant number of particles, 
there is no a priori reason to expect strong correlations 
between two specific ones. 

There are still many significant physics questions regard- 
ing jet energy deposition into the medium: Is there an 
appropiate theoretical framework to understand the jet 
energy deposition mechanism of a fraction / of the lost 
energy into collective modes in detail? Is sound the only 
collective mode compatible with the data? How can col- 
orful and colorless sound collective modes and their in- 
terplay be consistently included? How sensitive is the 
correlation pattern to details of the fireball evolution or 
the EOS? Can we successfully calculate the px spectra 
of near and far side simultaneously? At the moment, 
the case for the excitation and propagation of a sound 
mode through the hot and dense system appears to be 
rather strong and offers exciting opportunities to study 
the evolving system from a new angle, and we plan to ad- 
dress those questions in forthcoming publications along 
with a systematic investigation of some of our model as- 
sumptions. 
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